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ABSTRACT 


A numerical analysis of a heat transfer and combustion 
model for a porous medium within a circular cylinder is 
carried out. The basic mechanisms considered in the theory 
are: a carbon energy conservation equation, energy balance 
on the air, and an oxygen mass balance equation. Heat trans- 
fer mechanisms included in the model are conduction, convec- 
tion, and radiation. A heat generation term arising from 
combustion of the carbon is included in the model. Transport 
mechanisms for oxygen mass transfer are molecular diffusion 
and convective transport. The governing heat and mass trans- 
fer equations are solved numerically by the Galerkin formula- 
tion of the Finite Element Method. The results show the 
effects of permeability, porosity, geometry and initial condi- 


tion on the behavior Of “theme, stem. 


TABLE OF CONTENTS 


ii. TERR LE TSS) CONG AB SIE © | PIF a pe ala eZ 
eT. CYTES IRIE IE TON PIO IE ey 129810) 5) ELH i Oe le a td Zak 
TL BIOS OTE TE GLO JOSS MELD NS MSTA EL) SS aI el a tae a ar a 28 
OW PS eUoo ON WmeclemhUSlONS ===-—==-==---<<---------- 3: 
eo ee lower ene ollY AND PERMEABILITY --------- eel 
Ee ee Clo wOr wee EA SSS a = Se ee 44 
Sees, Cf Cou DAkY ~CONDITIONS -—-==-—---------- 49 
Dre ee eno Ohmi ela > = = a SJ 
Pec elem olON> Ane eke cCOMMENDATIONS: ====<-+=-------- 52 
Beene eee OU DARY CONDITIONS ==---==--==--=----------- 55 
en ® eee fs eee ee SS SS = = = - — = = = --------- 58 
fo LTE iby ID IE SAD ALY AUG BL JANG ta 60 


Ls 


ivi. 


iE ies 


LIST OF TABLES 


POROSITY AND PERMEABILITY -~--------------- 


POROSITY AND] PERMEABIEEIy FOR FIGURE 427 


POROSITY AND PERMEABILITY FOR FIGURE 4.8 


am aw em @ ow ow eo o= a= 


ioe OF FIGURES 


Geometric Model of Porous Medium ----------------- 
Grid Convergence for the Finite Element Method --- 


Carbon Temperature Profile with d/D = 1.0 and 
Constant Porosity (p = 0.476) -------------------- 


Oxygen Concentration Profile with d/D = 1.0 and 
Constant Porosity (p = 0.476) -------------------- 


Carbon Temperature Profile with d/D = 0.75 and 
Constant Porosity (p = 0.779) -------------------- 


Oxygen Concentration Profile with d/D = 0.75 
and Constant Porosity (p = 0.779) ---------------- 


Carbon Temperature and Oxygen Concentration vs. 
Time for Varying d/D ratioOS ------99- oreo nr nnn 


Carbon Temperature vs. Z2/z for Varying Time ------ 


Carbon Temperature and Oxygen Concentration vs. 
Time for Varying Permeabilities ------------------ 


Carbon Temperature and Oxygen Concentration vs. 
Time for Varying Porosities ---------------------- 


Carbon Temperature Profile with Varying Heat 
FLluX -occcc--- --- - - - - 


Oxygen Concentration Profile with Varying 
Heat Flux --------------- ---- -- - 


Carbon Temperature and Oxygen Concentration vs. 
ie on etic Nom eowmainy Comal t lon ——————————— 


Carbon Temperature and Oxygen Concentration vs. 
Time for Varying the Axial Length of the 
Senn eeu ee 


DD 


Qe We DP 


UO 


G 


AH 


Lh 


M 


LLST Ones MBOLS 


Arrhenius coefficient, wetted surface area 
Stiffness matrix 

Mass matrix 

Specific heat at constant pressure 
Unit cell dimension 

Mass diffusion coefficient 
Particle diameter 

Activation energy 

Specific internal energy 

Fraction 

Excitation vector 

Sto1chiometric racie@ 

Pseudo mass velocity 

Gravitational constant 

Heat of combustion 

Convection heat transfer coefficient 
Constante 

Thermal Conductivity 

Thickness of porous medium 

Element length 

Molecular weight 

Specific permeability 


Reaction order 


Pr 


ce. 


n-l 


V 


Z 


Greek 


Global Basis Function 
Eto S Seicgis 

Prandtl number 
Porosity 

Filter velocity 

Heat generation or loss 
Heat flux 

Reaction rate 

Rate per unit area 
Gas constant 

Reynolds number 
Radial coordinate 
Residual function 
Femperature 

Absolute temperature 
Time 

Previous time step 
Internal energy 

Pore velocity, radial. 
Pore velocity, axial 
Void volume, molecular volume 
Axial coordinate 


Specific internal area 


Symbols 


Thermal diffusivity 


Pore diameter 


6 Kronecker delta function 

E Thermal emissivity 

8 Solution coefficient 

. Field operator 

i Dynamic viscosity 

E Local element coordinate 

fe) Mass density 

O Stefan-Boltzman constant 

T Tortuosity, stress 

w) Oxygen concentration 

W Particle shape factor, approximate solution 

SUDS@GIOES 

a Air 

e Carbon 

CO Carbon monoxide 

CO, Carbon d10Oxide 

e Effective 

fm Film 

g Heat generation 

1 At the current time or step 

any ike pejihie ko; al 

O Cylinder dimension (i.e., r. 1S cylinder radius, 
Z 6 is cylinder length) : 

O5 Oxygen 

1 At COnstalweerecsure 

18 Radiation 


10 


Sst 


Solid 
S) Eeligiealigie) 
Universal 


Ambient conditions 


ale 


eo ObUCTLLON 


Combustion of porous media has far-reaching practical 
consequences. As there are large reserves of available eGau 
in many places, coal combustion iS a potential source of 
increased energy production if some of the pollution problems 
can be resolved. The combustion of coal is also useful for 
the heating of homes and buildings. On the other hand, unde- 
Sirable combustion of matter in the form of fires causes 
great losses in lives and property destruction every year. 

In order to gain the most advantage from its beneficial 
uses, and to minimize its destructive power, a clear under- 
standing of combustion behavior is necessary. Among the 
various engineering disciplines, combustion is one of the 
least understood. The combustion problem isespecially diffi- 
cult because of its multi-disciplinary nature, involving the 
areas of fluid mechanics, heat transfer and chemical kinetics. 
Chemical kinetics of combustion is a most complex subject in 
itself, and the difficulty is compounded by the additional 
complications introduced by the other disciplinary areas. 

Because of its importance, however, the area of combustion 
has been the subject of much investigation, both analytically 
and experimentally. In a 1973 survey paper, Emmons [Ref. 1] 
points out that the three major obstacles in the study of fires 


are turbulence, the immensely intricate chemistry, and the lack 
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of knowledge of the physical parameters dependence on tempera- 
ture and composition. Emmons discusses the essential role of 
heat and mass transfer in fires, and points out that the con- 
servation laws of mass, energy, and momentum must be supple- 
mented by many other facts and laws accounting for chemistry, 
feoeana liguid diffusion, mutual diffusion of gases, radiation 
effects, and the variation of thermodynamic, chemical, kinetic, 
transport, radiative and mechanical properties over the fire 
temperature range. Neverless Emmons assets "Now is the time 
to take a new look, to detect and analyze the simpler parts 

of a fire, and thus open up the future of rational approaches." 

Emmons points out that inclusion of the complete chemistry 
in a model would be overwhelming even if it were Known, and 
that only that chemistry which is essential to the problem 
be included in the formulation. A brief summary of several 
"simpler" problems provides some insight into the overall 
problem. 

In a 1936 paper, Parker and Hottel [Ref. 2], presented the 
results of an experimental study on the combustion of carbon. 
This eee cntall investigation considered cylindrical carbon 
specimens combusting in air. They determined the distribution 


of N COx CO,, and O, in the gaseous "film" surrounding the 


e 


carbon at varying ambient temperatures and air velocities. 
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Resulting combustion temperatures varied between 1000 K and 
IPOs 2526 Sto be 25°C, 1340°F to 2240°F). Results show that 


the partial pressure of CO decreases with distance fror the 
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carbon surface, while the partial pressure of 0. increases with 
distance from the carbon surface. At a fixed temperature, the 
partial pressure of 0. at the carbon surface increased with 
increasing air velocity; while at fixed velocity the pameuam 
pressure of oxygen decreased with increasing temperature. 
Additionally, a separate experiment measured combustion rate 
as a function of temperature and air velocity. Their results 
show that at constant velocity, the reaction rate increases 
with increasing temperature; and the reaction rate increases 
with increasing velocity when the temperature is held constant. 
The experiments demonstrated that there is little CO at the 
carbon surface at the higher temperatures. An equation for 
the reaction rate of carbon brush is presented. The denominator 
of this equation contains two additive terms, namely, the 
diffusional resistance Rg and the chemical resistance Ro: 
With increasing temperature, Ro decreases and the combustion 1s 
adiffusion controlled; when the velocity increases, Rs decreases 
and the combustion is kinetically controlled. 

In a 1951 paper, Arthur [Ref. 3] presents the results of 
an experimental investigation of the combustion of two carbons 
of widely different reactivities in the temperature range of 
735-1175 K (460-900°C, 860-1650°F). The study Shows thatwiwe 
reactions occur; one leads to CO. formation and the other to 
CO formation. It was shown that 1) the rate of carbon gasi- 


fication decreases but slightly with carbon consumption, 2) 


the combustion temperatures increase with air flow rate, 
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ey) the co/co., ratio increases with temperature, with cO/CO., 
ene aes Kk  (oos CC, 1533°R) and 11.1 at 1168 K (895°C, 
1645°F), 4) the gasification rate is independent of air flow 

rate at low temperatures and increases rapidly at higher 
temperatures (the transition occurred at 595 K for carbon A, 

ema 293 K for carbon B), and 5) the CO/CO. ratio was indepen- 
dent of the initial partial pressure of oxygen. The most 
Significant result is that in the 673 K to 1173 K temperature 
range, the CO/CO, ratio increases exponentially with temperature 
meom 13 at 673 K to 12% at 1173 K. 

In a 1956 investigation, Koizumi [Ref. 4] conducted an 
experimental and analytical study of combustion of solid fuels 
in fixed beds. He found 1) that the oxidation zone is only 
one or two partical diameters deep from the burning front, 

2) that the burning rate 1s controlled by the rate of gas 
diffusion in the boundary layer (i.e., the chemical reation 
occurs instantaneously, and 3) that the CO to co, combustion 
occurs in the boundary layer. 

In 1961, Green and Perry [Ref. 5] presented the results of 
an analytical investigation of heat transfer with a flowing 
fluid through a porous medium. They formulated a one-dimensional 
transient model including conduction and convection. The 
fluid flow rate was assumed constant, and the thermophysical 
properties were assumed to be independent of temperature. Two 


coupled p.d.e. were obtained from energy balances on the solid 


(porous media), and the flowing fluid. Boundary conditions 
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were taken as equal fluid and solid temperatures at X = 0, and 
X > », The initial condition was taken as equal fluid and 
solid temperatures for all X. The p.d.e. were solved by the 
finite-difference method. They found that heat conduction 

and convection are both important for internal Reynolds num- 
bers less than unity; and that the fluid and solid tempera- 
tures are essentially the same when the dimensionless parameter 


(nz/k,)*/*+ (k,/P_ u) 1S greater than 0.342. 


He 

In a 1974 paper, Anderson and Zienkiewicz [Ref. 6] pre- 
sented an analytical investigation of a reactive solid with 
zero order kinetics. Two problems were considered: the 
Steady state problem, and the transient problem. In the steady 
state problem, the Frank-Kamenetski parameter is determined. 
This parameter indicates whether combustion will occur or not. 
In the transient case, the body is immersed inathermal field 
whose temperature is greater than the ignition temperature of 
the body. Here the time to ignition (1i.e., the induction time) 
is determined, as well as transient solutions T(x,y,t) of 
two-dimensional bodies. The problem is solved by the Galerkin 
finite element method. Although the formulation permits varia- 
ble conductivity and snecific heats; the authors assume these 
parameters constant in their illustrative examples. A single 
nonlinear p.d.e. obtained by an energy balance, includes 
conduction and convection heat transfer models. 

In a 1976 paper, Kim and Chung [Ref. 7] presented the 


results of an analytical investigation of conditions leading 
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to a state of self-sustaining combustion in a porous medium. 
The one-dimensional transient model includes combustion and 
conduction. Two p.d.e. were obtained from an energy and 
mass balance of the porous medium. Thermophysical property 
dependence on temperature waS not accounted for. Laplace's 
asymptotic method was used to obtain a solution. Solid gas 
reactions are considered to occur volumically throughout the 
porous medium. A constant rate of oxidant fuel is assumed at 
X = 0. Zero gradients of temperature and oxidant concentra- 
tion at X + © were taken aS boundary conditions. The boundary 
conditions at X = 0 are -AOT/OX = q and -DdC/d9X = m. Their 
results for a semi-infinite slab show that ignition is strongly 
dependent upon temperature and reactant concentration. The 
Lewis number was taken as unity. They also show that as 
either the activation energy decreases or the initial tempera- 
ture of the porous medium increases, the ignition occurs faster. 
In 1976, Sawyer and Shuck [Ref. 8] presented a transient, 
one-dimensional model for underground coal gasification of a 
coal bed with vertical micro fissures along which combustion 
occurs. Mass and energy balance for the solid and gas phases, 
and mass balance equation for each Species (oxygen, carbon 
monoxide, and carhon dioxide) leads to 7 coupled, nonlinear 
p.d.e. The heat transfer mechanisms include conduction in the 
solid, convection in the gas and convection between the solid 
and gas. All properties except permeability and reaction rates 


were assumed constant. The 7 p.d.e. are supplemented by an 
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equation of state for the gas, and a Darcy flow equation in 
the fissures. The system is solved by the finite-difference 
method. 

Sahota and Pagni [Ref. 9] present the results of an analy- 
tical investigation of heat and mass transfer in porous media 
subject to fires in a 1978 paper. The one-dimensional tran- 
Silent model includes conduction and convection. Energy and 
mass conservation equations, together with species conserva- 
tion equations, Darcy's law and state equations lead to 8 
equations which are solved by the finite difference method. 

A portion of the domain has tae: moisture. The energy 
balance includes a heat sink term to account for the evanora- 
tion of water. Thermophysical properties are treated as con- 
stants. When pressures are not required, a simplified technique 
for determining temperature has been developed. The method 
neglects heat transfer by convection and mass diffusion. This 
analysis is valid when very small air vapor mixture velocity 

1s present, and for very small Lewis numbers (1.e., when 

era << and ORT ETS oe Le 

In a 1980 paper, Saatdjian and Caltergirone [Ref. 10] pre- 
sented a mathematical model for a plane porous matrix bounded 
by two impermeable surfaces. The porous matrix 1S composed 
of a mixture of nonreacting fiber glass and a resin which 
decomposes exothermically into gaseous products as the matrix 
temperature increases. The two-dimensional transient model 


includes conduction, convection and combustion, but does meet 
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account for temperature dependent properties. The fluid is 
modelled as a perfect gas, and flow through the porous matrix 
is governed by Darcy's law. It is assumed that the matrix 
properties remain unchanged during the reaction (i.e., the 
fiber fraction/resin remains greater than 0.90), and the 
gaseous products and the fluid which initially filled the 


enclosed region are the same. The density of the sSaturating 


gas is allowed to vary. Results show temperature and density 
as functions of time. The problem counvles the effects of 
two phenomena, convective movement and chemical reaction. As 


the reaction is depleted, natural convection alone subsists. 
In a 1981 paper, Chan and Banejess [Ref. 11] investigated 
transient three-dimensional natural convection in porous media. 
The model includes a continuity equation, two momentum equa- 
tions and two energy balance equations. The finite difference 
method solves the five p.d.e.'s for the solid and fluid tem- 
peratures, the velocity components in the X and Y directions, 
and the pressure. The method uses upwind differencing for 
the advective term. This model allows for unequal solid and 
fluid temperatures. The authors state that the solid and 
fluid temperatures are equal ("equal diffusivity" model) when 
the filtration velocity is not too high and both phases are 
well dispersed. Fictitious conductivities take into account 
the dispersed structure of the solid matrix and the effect of 
hydrodynamic dispersion in the fluid energy balance equation. 


The porous medium is modelled as an assembly of spheres. 
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The authors introduce a scheme to reduce diffusional trunca- 
tion errors. Analytical results show good agreement with 
experimental results for Rayleigh numbers between 36 and 800. 
Numerical results at higher Rayleigh numbers deviate from 
experimental results. The model predicts convective insta- 
bilities at high Rayleigh numbers. 

In a 1981 paper, Hickox [Ref. 12] presented an investiga- 


tion into convective flow, induced by a heat source, through 


a porous medium of low permeability. The two-dimensional 
transient model includes conduction and convection. Four 
equations are obtained: an equation of continuity, Darcy's 


law for fluid flow, an energy balance on the porous medium 

and a state equation. The problem arises, for example, when 

a container of radioactive waste material 1s deposited ina 
sea bed below the sedimentary layer. The BousSinesq approxi- 
mation that density changes are due only to huoyancy is used. 
Permeability, viscosity, thermal conductance, thermal capacity 
and the thermal expansion coefficient are treated as constants. 
Numerical solutions based on similarity transformations were 
obtained for steady state point and line sources, as well as 

a steady state, constant temperature sphere. The results 

are predicted to be accurate for Reynolds numbers less than 
unity. In the case of the subseabed porous medium, a Reynolds 
number of 10 1s anticipated. The analysis is valid only for 


distances a few cannisters away from the disposed cannister. 
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ie CE SCRIPTION OF THE PROBLEM 


This investigation considers the combustion of a porous 
graphite medium. The porous medium is constructed from 
graphite particulate matter, either in the form of spherical 
particles or cylindrical fibers, imbedded in a cylindrical 
container. We assume the porous medium has a regular periodic 
structure with interconnecting pores which permit the flow of 
air through the medium. Figure 2.1 shows a geometric model 
of a regular periodic structure. 

Initially the porous medium is at ambient conditions. 

It 1s then subjected to a heat source at the bottom of the 
vertical container for a specified period of time. Upon 
termination of the heat flux, the system is at a greater 
temperature (and lower oxygen concentration) than the ambient 
temperature. Thereafter the system temperature will either 
decrease to ambient temperature, or it will continue to burn 
at higher temperature. The former process leads to extinction, 
the latter process is combustion. Which of these processes 
occurs depends on the particular attributes of the system such 
as its porosity, permeability, thickness and other distinguish- 
ing properties. This investigation seeks to determine how 
combustion of the porous graphite medium is related to the 
properties of the medium. 

The air which flows through the porous medium has two 


effects: it supplies oxygen for combustion which produces 
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Figure 2.1 Geometric Model of Porous Medium 


heat, and it provides for convective heat transfer through 
the medium. Combustion occurs if the heat generation of 
combustion is the dominant effect. If the heat transfer is 
meme aominant effect, then extinction will occur. Which of 
these effects dominates depends upon the physical system 
under consideration, as well as the boundary and initial 
@enditions. 

The mathematical model is formulated in the following way. 
We assume the air flow through the porous medium is governed 
by Darcy's law. Previous investigation has shown that the 
Darcy law model is an excellent approximation when the Reynolds 
number is less than 1. This occurs for the flow through the 
porous media covered in this investigation. In any case, that 
Reynolds number is less than one must be verified from the 
results of the analysis. The use of Darcy's law obviates 
the conservation of momentum equations. This results ina 
very Significant simplification of the mathematical model, 
with a corresponding reduction of computational effort. 

Energy balance on the porous solid, and on the air flowing 
through the porous solid vield two heat transfer equations. 
These heat transfer equations include conduction, and convec- 
tion mechanisms. In addition the porous solid energy equation 
includes a radiation term; and the air energy baJance equation 
includes a term to account for the energy transport by flow. 
The porous solid energy balance equation also contains a 


heat generation term to account for the combustion of the 
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porous graphite. This term is obtained as the product of the 
reaction rate and enthalpy of formation. Both energy balance 
equations have a term reflecting the change in internal energy. 
We note, in passing, that the independent variables are the 
Spatial cylindrical coordinates r and z, and the time coordinate 
t. The dependent variables are the temperatures of the porous 
graphite Ty: and the air flowing through the medium me In 
addition the heat generation term in the porous solid energy 
balance equation introduces another dependent variable, the 
oxygen concentration, into the equation. 

The third field equation is obtained by a mass balance 
of oxygen. This oxygen transport equation includes a molecu- 
lar diffusion term, a convective transport term and a term 
to account for oxygen consumption due to combustion. fThis 
last term is obtained as the product of the reaction rate 
and the inverse of the stoichiometric ratio of the reaction. 
Diffusion due to temperature and pressure gradients were 
assumed to be negligible. 

Finally a fourth field equation was obtained by combining 
Darcy's law for flow with the continuity equation. This 
equation introduces pressure as the fourth field variable. 

The above description 18S a brief summary and is not com- 
plete in and of itself. Each of the individual terms in 
each of the field equations must be modeled in their own right. 
For example, there is a model for the diffusion coefficient 


of oxygen into air; there 1s a model for the effective 
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Sencquctivity parameter for the porous solid, and so on. This 
investigation did not develop each of these models, and there- 
fore the interested reader can obtain details of all of the 
models in the original developments of Vatikiotis [Ref. 12] 
and Martinez [Ref. 13]. 

For the sake of completeness, a brief description of the 
combustion and chemical kinetic models follows. The essen- 
tial ideas of the combustion model were presented by N.N. 
Semenov [Ref. 14]. Semenov's model proposes that the combus- 
tion process occurs in two phases, the kinetic regime, and the 
diffusion regime. The initial combustion at lower temperatures 
is called the kinetic regime because there is sufficient oxy- 
gen present for combustion to proceed and the reaction is 
controlled by the actual kinetics of the surface reaction and 
does not depend on the diffusion rate. In the kinetic re- 
dime, the reaction rate increases, exponentially with tempera- 
ture. As the temperature increases, the reaction rate increases 
and the rate of osygen consumption increases; causing a lack 
of oxygen for further increase in the reaction. In this case, 
the lack of oxygen slows down the reaction. This regime of 
combustion, whereby the reaction rate is limited by the diffu- 
Sion of oxygen is known as the diffusion regime. A graph of 
reaction rate versus temperature results in the well-known 
S-curve. The lower portion of the S-curve defines the kinetic 
regime; the upper portion of the S-curve defines the diffusion 


regime. 
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The chemical Kinetics of combustion is much too complex 
to model all chemical reactions in detail. Instead, an 
approximate model was used whereby we assumed that the products 
of graphite combustion were only carbon monoxide and carbon 
dioxide. Empirical results by other investigations have 
proposed such models. In the present work, the model pro- 
posed by Arthur [Ref. 3] was taken. This model presents an 
expresSion which relates the ratio of CO and CO, to the 
temperature of the reaction. 

The expresSion used for the reaction rate itself, is of 
the Arrhenius type. However, in accordance with the remarks 
of Frank-Kamenetskii [Ref. 15], the order of the reaction was 
not set to unity. Frank-Kamenetskil proposed that a fractional 
order reaction between 1/3 and 2/3 is a better representation 
of the experimental data of Parker and Hottel [Ref. 2]. This 
parameter was left to the selection of the investigator. It 
should be noted that combustion analyses show that significant 
changes in the system behavior occur with small changes in the 
reaction order. 

The final mathematical model consists of four transient 
field equations in the four unknowns: graphite temperature, 
air temperature, oxygen concentration and pressure. These 
partial differential equations are both coupled and nonlinear. 
Together with boundary and initial conditions, these equations 
form a well-posed problem which cannot be solved exactly. An 


approximate solution was obtained by transforming the set of 
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4 p.d.e. into a system of ordinary differential equations 
uSing a Galerkin formulation of the finite element method. 


Details of the finite element formulation can be found in 


Peer. 13]. 
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III. NUMERICAL CONSIDERATIONS 


A Galerkin formulation of the finite element method was 
used to obtain solutions of the porous medium and air energy 
equations, the oxygen diffusion equation, and the continuity 
equation. 

A grid arrangement as shown in Figure 3.1, which has 7 


non-uniform divisions in the r—-direction and 14 non-uniioms 


divisions in the z-direction, was selected. The z-direction 
divisions were smaller near the Z/Z. = 0.0, where the carbon 
temperature and oxygen concentration change rapidly. The 


r-direction divisions were uniform in the middle, and non- 
uniform at the two ends a1 = 0.0 and r/x, = 1.0, where the 
Carbon temperature and oxygen concentration change slightly. 


The r and z grid lines are as follows: 


r/r,? O2G00, 0.125, 90.250, U-3500, Urt50, Weer 


OO 
Z/Z 5° 0. 000,—-0:..001,.. (0.0025, 92005, (050075, 20 7e 
0.015, 0.020, "0270500, 02100, UF 2500) Soe ten 


O70; eee 


The results obtained from this grid compared very well (within 


5% at all nodal points) with results obtained by a finer 
14 x14 grid, and therefore the 7 x14 grid was used for all 


computer analyses. 
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Figure 3.1. Grid Convergence for the Finite Element Method 


25 


The inputs for the integration program include (1) initial 
time step, (2) minimum time step, (3) maximum time step, and 
(4) property up-date time. The selection of these parameters 
affects the accuracy of the numerical solution. However, it 
is not possible to determine the error associated with any 
particular selection of these parameters. The selection of 
these parameters is by trial and error. If they're made 
very "small" the program runs very slowly and meaningful 
results could be very computationally expensive to obtain. 

If these parameters are selected too "large" the program may 
give poor resultS or it May not run at all. For this progeaie 
the initial time step is ios hr, the minimum time step is 


=o 
106, hr, the maximum time step 1s 0.25 hr, ana@ the up ldag. 


; ; -14 ae . 
time is 10 hr. If the initial time step and the up-date 
time is larger than a certain time interval, the program 


overflows or underflows in subprogram NUITSL and LDASUB. 
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ie SISCUSsONeAND CONGLUSIONS 


A heat transfer and combustion model for a porous medium 
within a cylinder is considered. The assumption of the model 
is that of the idealized porous medium described in Chapter 
III. The results of this investigation are intended to des- 
cribe the behavior of a porous carbon medium with a heat 
generation term. The behavior of the temperature and oxygen 
concentration of the system depends upon the initial conditions, 
the boundary conditions, the geometry of the system, the 
porosity, the permeability, and the heat flux. A number of 
analyses were performed to determine how some of these param- 
eters affect the system behavior. The boundary Donations 
for the investigations of Sections A, B, C, and D are pre- 
sented in the Appendix. The porous media wasS initially ata 
constant temperature of 100 degrees Fahrenheit, a constant 
meoyocen concentration of 0.0172 Lbm/£t? (ambient conditions), 


and a constant pressure of 2118.7 iy yaya Grae] OS1). 


Peres PheCTS OF POROSITY AND PERMEABILITY 

Section IV.A.1 presents the results of varying d/D ratios 
Miable TL, Bigures 4.1—-4.6) ..section IV.A.2 changes D and d, 
keeping the porosity constant, with varying permeability 
Stole 11, Figure.4.7). In Section, IV.A.3, the d/D ratio is 
varied so that permeability 1S constant and porosity is 


changing (Table III, Figure 4.8). 


ok 


1. Etiects*ofeeme ed Dekat ro 

In this section the effects of porosity and permeability 
on system behavior were investigated. Varying the carbon 
diameter, all other parameters fixed, results in changes to 
both the porosity and the permeability. Table I shows the 
three values selected for the carbon diameter and the 
resulting values of porosity and permeability. The results 
of this investigation are shown in Figures 4.1-4.6. In all 
investigations in this section a heat flux of 500 Bue eaeae 
was applied for 30 minutes and then turned off. Figures 4.1 
and 4.2 show the temperature and oxygen concentration profile 
with d/D = 1.0. Figures 4.3 and 4.4 show the temperature and 
oxygen concentration with d/D = 0.75. Figure 4.5 compares 
the carbon temperature and oxygen concentration as a function 
of time, at the posttioen L/L = Oe: and Z/Z., = 0007 scene the 
three sets of d/D ratios. Figure 4.6 shows how the carbon 
temperature varies with the z-coordinate axis and time, at 
r/X = 0m oe 

As shown in Figures 4.1-4.4, the d/D ratio affects the 
carbon temperature and oxygen concentration. These figures 
show that these variables are essentially indevendent of the 
radial direction, as they should be for the boundary condi- 
tions of these problems. Figures 4.1 and 4.2 show that at 12.2 
minutes and Z/Z = 0.0, the carbon temperature is 480 degrees 
Fahrenheit, and the oxygen concentration has not changed 


because the temperature ts too low. However, we note that the 


BY 


TABLE I 


POROSITY AND PERMEABILITY 
Pe \\, ‘ 


a/D = 10 a/D = 02873 qd/D = 0.75 


shape 


mit Cell 
Dimension 
Cit) O60 04 17 O [0004 ©0004 17 


Particle 
Diameter 
Crt ) O200c 0417 OC UO0s65 020003 13 


Porosity 0.476 0.649 G..779 


Pore 
Diameter 4 ‘ 


(£t) 30395 LON 5,005 <10- ene. Wap 


Specific 

Internal 

Area 

ee rt) 2) Taye. S) fate tse 10) Bele. 9 


Permea- 
eon Ja og 


(£t2) Saye eee 9 


2 


Pears. 10 8.959 x10. 
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1000. BUT/FT2-HA. 


CARBON TEMPERATURE SURFACE 
PROBLEM TIME i$ 6.10 MIN. 
HEAT RUX AT Z/ZO @ 

RO= 100 FT. 

ZO @ 1.00 FT. 


TMAX = 620. DEG. F. 


HEAT RLUX AT Z/ZO @ 600. BUT/FT2-HA. 


RO= 100 FT. 
ZO@= 1.00 FT. 


CARBON TEMPERATURE SUAFACE 
PROBLEM TIME IS 5.20 MIN. 
TMAX = 363. DEG F. 
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CARBON TEMPERATURE SURFACE 


CARBON TEMPERATURE SURFACE 


PROBLEM TIME I$ 420 MIN. 


PROBLEM TIME I$ 6.00 MIN. 
HEAT FLUX AT 2/ZO = 
RO= 100 FT. 

ZO = 100 FT. 


00. BUT/FT2-HA. 


HEAT FLUX AT Z/Z0 @ 2000. BUT/FT2-HA. 


RO= 100 FT. 
ZO= 100 FT. 


TMAX @ 


1001. DEG. F, 


TMAX = 886 DEG. F. 





VATDAIQAY AWBNLWWAAWAL 


(p = 0.476) 


Carbon Temperature Profile with d/D 


and constant porosity 


Figure 4 ou 
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02 CONCG.\® LOEFA LBM. IETS) 


O02 CONCENTRATION SURFACE 

PROBLEM TIME IS 620 MIN. . 

HEAT FLUX AT Z2/ZO = 600. BUT/FT2-HR. 
RO = 1.00 FT. 

ZO = 1.00 FT. 

O2 CONC, MIN = 172. 





O2 CONCENTRATION SURFACE 

PROBLEM TIME IS 5.00 MIN. 

HEAT FLUX AT Z/ZO = 1600. BUT/FT2-HR. 
RO = 100 FT. 

Zo 1.00 Fi: 

O2 CONC, MIN = 119. 
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O02 CONC. AOEF4 LEW. IETS) 


O2 CONCENTRATION SURFACE 

PROBLEM TIME !IS_ 5.10 MIN. 

HEAT FLUX AT Z/ZO = 1000. BUT/FT2-HA. 
RO = 100 FT. 

ZO = 1.00 FT. 

O2 CONC, MIN = 172. 





O2 CONCENTRATION SURFACE 

PROBLEM TIME IS 4.20 MIN. 

HEAT FLUX AT Z/ZO = 2000. BUT/FT2-HR 
nO =reo rl. 

ZO = 1.00 FT. 

O2 CONC, MIN = 0. 





Figure 4.2 Oxygen Concentration Profile with d/D = 1.0 


and COnstamt COLrosi ry 
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(p = 0.476) 


HEAT FLUX AT Z/ZO = 600. BUT/FT2-HA 


RO = 1.00 FT. 
ZO = 100 FT. 


CARBON TEMPERATURE SURFACE 
PROBLEM TE i$ 220 MIN. 


LITTLE O/ BG O= 0.76 


TMAX @ 


CARBON TEMPERATURE SURFACE 
PROBLEM TIME 18 0.20 MIN. 

HEAT FLUX AT Z/ZO @ 500. BUT/FT2-HA. 
RO= 100 FT. 

ZO = 100 FT. 

UTTLE OY BG DO = 0.75 


627. DEG. F. 


TMAX = 278. DEG F. 


CARBON TEMPERATURE SURFACE 
PROBLEM TIME IS 12.20 MIN. 
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100 FT. 
ZO= 100 FT. 


HEAT FLUX AT Z/Z0 @ 600. BUT/FT2-HR. 
LITTLE DY BAO = 0.76 
TMAX = 1202 DEG F. 


RO = 


600. BUT/FT2-HA. 
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Figure 4.3 


Carbon Temperature Promien witha, 


and constant porosity 


Ory 7 9) 
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O02 CONC.\S VDE A LBW 


O2 CONCENTRATION SURFACE O2 CONCENTRATION SURFACE 


PROBLEM TIME IS 0.20 MIN. PROBLEM TIME IS 220 MIN. 

HEAT FLUX AT Z/ZO0 = 600. BUT/FT2-HR HEAT FLUX AT Z/ZO = 600. BIZT/FT2-HR. 
nO = FOC Fil. RO = 100 FT. 

ZO™= 1.00 FT. ZO™ 100 FT. 

UTTLE D/ BIGD = 0.75 LITTLE D/ BIGD = 0.76 

O2 CONC, MIN = 172. O02 CONC, MIN = 172. 
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O2 CONCENTRATION SURFACE Q2 CONCENTRATION SURFACE 

PROBLEM TIME IS 5.00 MIN. PROBLEM TIME IS 12.20 MIN. 

HEAT FLUX AT Z/Z0 = 600. BUT/FT2-HR. HEAT FLUX AT Z/Z0 = 600. BUT/FT2-HA. 
nO = POO oT. RO = 100 FT. 

Zo = 0G FT. ZO = TeO ii, 

Bice ean = 0.76 Piiiee ©) sia b= 0.76 

QO2 CONC, MIN = 167. O2 CONC, MIN= 0. 


O02 CONC. VQEX4A LBW .PET3) 





Figure 4.4 Or seme COnecmemertOnmurro;1 le with d/D = 0.75 
BiGweOM Scant OOTOstty 4p = 0.779) 


a7 


Solely peourate Ae Los 
ETC me) Weakxg PUP ainieriedusj, uoqie9 G* p ealns’ty 


440 XN19 Lv 


soe oe 4IO XNA LWGH 
$ Oia \ (saat) Jet 
ce Ge ce ct i) 
eae aliregat aie ere o 
i\ | 
r ’ ‘ 8 3 
i i ' 4 
: : \ 
yee \ ‘ % 
ee \ S10 = oa/a\ 
r : 
ae ‘ | 8 
-—" % ; i 
$180 = oa/at\ \ » e C280 = Qa/d1 - 
‘\ 4 eg 
ae . - i at i. 
ie ca/oi ae ‘. : & 1 


3,001 = “dN3L WUIN 
00 = 022 
SO = ON/H * NOUWISOd 
YH-ZLI/NLA OOS = XN LV3H 


ee ee 





(4 JayNLva3G NIL 


Bye: 


auty sutAdep doj 


010 


4 00l = ‘dIV31 TWIN 
SO = Ou/e * NOUISOd 
YH-ZL4/NB OOS = XML LV3H 


6/90 = OFYM 


ine: 


/ 


L 


oS GN 


Ww) 


A oanjededuay, 


620 = OF/M 





O10 90°0 


uo Gu Lege) 





9° Ff 


aInNst J 





Se 


temperature gradient is very large over a small region in 
the vicinity of the imposed heat flux. The region of tempera- 
ture activity does not extend beyond Z/Z = 0.1. Figures 
4.3 and 4.4, with d/D = 0.75, show the carbon temperature and 
oxygen concentration. At 2/2, = 0.0 and 12.2 minutes, me 
carbon temperature is 1200 degrees Fahrenheit and the oxygen 
concentration becomes almost zero. The results show that 
combustion activity increases with decreasing d/D. This can 
be attrioutea™ to ain imeredse intimal owe eloeie 

Figure 4.5 shows the effect of d/D ratios on combustion 
aCtivalty at r/r = 0.5 and 2/2. = 0.0. As shown in Figure 4.5, 
with @/D = 1.0, combustion activity is the least and the oxygen 
remains almost constant. With d/D = 0.875, combustion activity 
increases so that the oxygen is depleted at Z/Z = 0.0 in 
about 30 minutes. With da/D = 0.75, combustion activity becomes 
very strong resulting in a depletion of oxygen at Z/ 2. = 0.0 
in about 10 minutes. The reason being that as the d/D ratio 
decreases, there 1S an increase in void space, which allows 
more oxygen to be supplied to the combustion region. Also, 
with a smaller d/D, the carbon mass per unit cell is smaller 
and therefore the heat flux per carbon mass increases. This 
results in an increase of the carbon temperature. In the case 
of ad/D = 0.75, there was no opportunity to monitor behavior 
beyond 10 minutes because the program terminated due to numeri- 
cal difficulties. For the other two cases of @/D = 0.875 


and 1.0, the heat flux was removed after 30 minutes to observe 
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whether the system would continue to combustion Orgo to extinc- 
tion. The carbon temperatures dropped sharply and the porous 
carbon medium returned to its initial temperature while the 
oxygen concentration returned to its original concentration 
somewhat less sharply. 

Figure 4.6 illustrates the behavior of carbon tempera- 
ture versus 2/2. at r/T = 0.5 as time 1S varied for each of 
the d/D ratios. These figures show that the carbon tempera- 
ture 1S greatest at Z/ 2. = 0.0, and the carbon temperature 
approaches the constant ambient temperature as 2/2 increases. 
Thus in all cases we observe that combustion occurs in a local 
region near the base of the system where oxygen is entering 
into the medium. 

°2. Effects of Permeability 

Figure 4.7 shows the carbon temperature and oxygen 
concentration as a function of time, at r/X =a. 5 and 2/2. 
= 0.0, for various permeabilities with porosity fixed. The 
boundary conditions are the same as Section V.A.1 except that 
the heat flux was changed from 500 Btu/f£t--hr to 1000 Bey Be 
hr. The description of the porous medium used in this inves- 
tigation is shown in Table II. Porosities and d/D ratios are 
the same for cases A, B, C, and D. Case A has the smallest 
unit cell (D = 0.0002085 ft), the smallest particle size, and 
the smallest permeability (6.993 eG eee case D has the 
largest unit cell (D = 0.000834 ft), the largest particle size, 


and the largest permeability (1.119 lon BE). Unfortunately 
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TABLES. 


POROS Ii, = >eeERMEABILETY FOR FIGURE 4.7 


Unit Cell Particle Ratio Porosity Permeability 
Dimension Diameter (ft-) 

(Et) Gegey: 
.0002085 .0001824 .875 0.6492 5 Cee sany oo 
.000417 .0003649 875 0.6492 De am 
.0006255 .0005473 875 0.6492 6.29 x107° 
000834 0007298 875 0.6492 at, oe 
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it was not possible to obtain analysis beyond t = 5 minutes 
because the program has difficulty running when the oxygen 
concentration gets very small. The results show that as 
permeability increases, the carbon temperature slightly 
increases, while the oxygen concentration increases. This 
shows that an increase in permeability provides additional 
oxygen for a small enhancement of combustion. 

3... Effects Of Pesesicy, 

In this section the effects of porosity on system 
pehausere are investigated. Figure 4.8 shows the carbon 
temperature and the oxygen concentration, at r/T = 0.5 and 
2/2. = 0.0, as a function of time. The description of the 
porous medium used in this investigation 1s shown in Table 
III. Here porosity is the same for cases A and B, and cases 
C and D, while permeability is the same for cases A and C, and 
cases B and D. Fiaure 4.8 shows that for the same porosity 
(cases A and B, and cases C and D) carbon temperature and 
oxygen concentration are the same. For the same permeability, 
the carbon temperature and oxygen concentration have different. 
values. It appears from these results that the carbon tempera- 
ture and oxygen concentration are strongly affected by porosity 
and are not affected by permeability. This result is diffi- 


cult to explain and may indicate an error in the program. 


rae BPRPECTS OF (Nba hia 
Figures 4.9 and 4.10 show the changes in carbon tempera- 


ture and oxygen concentration as the heat flux is varied. 
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TABLE III 


POROSITY AND PERMEABILITY FOR FIGURE 4 .¢ 


UniteCert Particle The d/D Porosity Permeability 
Dimension Diameter Ratio (Ft 

(ft) (ft) 
000866 000866 1.0 0.4764 >. 797 sane 
001733 001733 1.0 0.4764 1.119) Sia 
.000417 .000365 0.875 0.6492 2.797 x10” 
000834 000730 0.875 0.6492 1,119 ~ toee 
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1000 BUT/FT2-HR. 


CARBON TEMPERATURE SURFACE 
PROBLEM TIME !$ 6.10 MIN 
HEAT FLUX AT Z/Z0 = 

RO= 100 FT. 
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CARBON TEMPERATURE SURFACE 
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02 CONCENTRATION SURFACE 

PROBLEM TIME IS 620 MIN. 

HEAT FLUX AT Z/ZO0 = 600. BUT/FT2-HR 
RO = 1.00 FT. 

ZO = 100 FT. 

O02 CONC, MIN = 172 


O2 CONCENTRATION SURFACE 

PROBLEM TIME IS 600 MIN. 

HEAT FLUX AT Z/ZO = 1500. BUT/FT2-HR. 
RO = 100 FT. 

ZO= 1.00 FT. 

O2 CONC, MIN = 119. 
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O02 CONCENTRATION SURFACE 

PROBLEM TIME IS_ 5.10 MIN. 

HEAT FLUX AT Z/ZO = 1000. BUT/FT2-+R. 
RO = 100 FT. 

ZO = 1.00 FT. 

Q2 CONC, MIN = 172. 





O2 CONCENTRATION SURFACE 

PROBLEM TIME IS 420 MIN. 

HEAT FLUX AT Z/ZO = 2000. BUT/FT2-HR 
RO = 100 FT. 

ZO = 1.00 FT. 

O2 CONC, MIN = 0. 
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Heat flux was varied 4 times (500, 1000, 1500, and 2000 
eeeset--hr) . 

The relationship between the carbon temperature and the 
heat flux, after exposing the system to about 5 minutes of 


heat flux, was observed as follows: 


Heat Flux Heated Time Carbon Temperature 
fptn/ftsenr) (minutes) (degrees F) 
500 O22 B10 
1000 Det 620 
L500 5 0) 885 
2000 4.2 Oa 


Heat fluxes beyond 2000 a ae cause the oxygen concen- 
tration to be depleted before 5 minutes, making the data 
unusable for any valid comparisons. 

The relationship between the oxygen concentration and 
heat flux, after exposing the system to approximately 5 


Minutes of heat flux, was observed as follows: 


Heat ae Heated Time Oxygen Concentration 
eu / Eon, ) (minutes) (1bm O,/ft ) 

500 3.0 192 ne 

1000 Sil Wo aes 

1500 5.0 mo) Vom 

2000 4.2 0 


ee PPePECLS OF BOUNDARY CONDITIONS 
Figure 4.11 shows the results of the temperature and 


oxygen concentration with the same initial conditions and 
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boundary conditions as before, with a single exception. In 
case B, the boundary condition for the oxygen concentration 
equation at Z/Z = 0.0 was changed from a Cauchy boundary 
condition to a constant oxygen concentration boundary condi- 
tion. As seen in Figure 4.11, een temperature and oxygen 
concentration depend on the boundary condition. The results 
Show that whereas the case B boundary condition provides a 
constant supply of oxygen at Z/Z = 0.0, the case A boundary 
condition leads to a decrease in oxygen, starting about 

t = 4.5 minutes, at that boundary. For both cases (A and B), 
the carbon temperature remains the same until the carbon 
temperature reaches 860 degrees Fahrenheit which occurs at 

t = 12.0 minutes. Thereafter, the carbon temperature for case 
B increases faster anon the carbon temperature for case A. 
This may be attributed to the additional oxygen at the boundary 


associated with case B. 


Pee EFFECTS OF GEOMETRY 

Figure 4.12 shows the carbon temperature and oxygen concen- 
fiat LOn sat Z/Z = 0.0 and E/E SO = 085, for three cylinders of 
different axial lengths. The heat flux for each case is 1000 
pen fe hw Varying the axial length of the cylinder, the 
carbon temperature and oxygen concentration yields different 
results. As shown in Figure 4.12, the carbon temperature in- 
creases with increasing axial length, while the oxygen concen- 


tration decreases with increasing axial legnth. These results 


ai 


ZIPUTTAD Sa FO Yabusy Terxy eya burAaeA 703 
OWT], “SA *‘ouo0g UussdsAxQ PUP dUN ,edJaduay UOGUeD 


(S14/NEI)ONOD NGOAXO 





43,001 = “dA3l “LN 

00 = O2/2 
SO = O4/Y -NOWSOd 
YH-714/NLG OOO! = XN LY3H 


ZL’ p eansty 





De 


are difficult to understand, and may indicate that the program 


ieenet perctorming correctly. 


E. CONCLUSIONS AND RECOMMENDATIONS 

This analysis shows that the behavior of combustion and 
heat transfer in a porous carbon medium is a complex process. 
The behavior of the system depends upon porosity, permeability, 
boundary conditions, initial conditions, and geometry. Some 
@mueene results of this investigation could be satisfactorily 
expected (the effects of the d/D ratios, and the boundary con- 
ditions), while others could not (the effects of permeability, 
and the effects of the geometry). The computer program developed 
by Lt. Martinez was used. This computer program was initially 
confirmed by Lt. Martinez by the model validation tests on 
cylinder problems with known solutions. However during this 
investigation, two difficulties with the computer program were 
encountered. First, the computer program terminates when a 
carbon temperature is higher than 1300°F and an oxygen concen- 
tration is less than nome inven This indicates that the 
integration algorithm was not able to perform satisfactorily 
when severe gradients occurred. Secondly, some results, such 
as the effects of permeability change and geometry change, 
are unexpected and infer that the program may have some errors 
edit . 

Improvements may be realized by extending the porous 


eelsoon medium mogel to include: 


jo) 


ee) 


development of an integration algorithm which is 
capable of performing satisfactorily over a larger 
range of temperatures; 


a detailed investigation of the program to determine 


whether there are errors associated with permeability 
and geometric effects. 
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APPENDIX 


BOUNDARY CONDITIONS 


. Boundary conditions employed were as follows for the 


Garpon, 
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Monally, for 0.5 concentration, 
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